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Abstract 

We study the correlations (and alignment as a particular case) existent 
between the fragments originated in a decaying process when the daughter 
particles interact. The interaction between the particles is modeled us- 
ing the potential of coupled oscillators, which can be treated analytically. 
This approach can be considered as a first step towards the characteri- 
zation of realistic interacting decaying systems, an archetypal process in 
physics. The results presented here also suggest the possibility of manip- 
ulating correlations using external fields, a technique that could be useful 
to provide sources of entangled massive particles. 
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1 Introduction 

In classical mechanics two particles originated in a decaying process maintain 
perfect correlations between their positions and momenta during the subsequent 
evolution of the system. In quantum mechanics the situation is not so favorable. 
The spread of the wavefunction and the uncertainties in the initial positions and 
momenta can cause in many cases a rapid lost of the initial correlations [1, 2, 3]. 
In these cases the correlations are not useful for practical purposes. References 
[2, 3] studied free non-relativistic decaying systems. Real systems are not in 
general free, but present various types of interactions. We pose in this paper 
the question if, when these interactions are taken into account, the picture of 
the process and the behavior of the correlations, are modified. The study of real 
interactions is in general very complicated, and many times we must resort to 
perturbative methods. We propose here an alternative approach. We consider 
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a simple model of interaction that can be solved analytically by assuming that 
the two particles interact as coupled oscillators. In all the areas of physics, but 
in particular in quantum mechanics, many problems that are not analytically 
treatable are reduced to harmonic oscillators. As we are dealing with two par- 
ticles, coupled oscillators could be a good candidate for our problem. On the 
other hand, coupled oscillators have been used in several branches of physics to 
simulate more complex systems [4, 5, 6]. We remark that we are modeling the 
evolution taking place after the decaying process, not the transformation of the 
mother into two daughter particles, a non-unitary process incompatible with 
the unitary evolution of coupled oscillators. The model of interaction not only 
holds for the direct interaction between the two particles (the electromagnetic 
interaction in the case of charged particles) but also for external fields, which 
could be introduced with the aim of driving the evolution of the system (see, for 
instance, Ref. [7] where external fields are used to prepare non-spreading wave 
packets). 

In addition to the more complete characterization of decaying systems (one 
prototype example in physics) provided by the inclusion of interactions, the 
analysis below allows for a better understanding of how correlations could be 
manipulated. In connection with the last point and from a more practical point 
of view, it must be remarked that decaying systems with good correlations could 
provide a source of entangled massive particles, well suited for two-particle in- 
terference experiments. Up to now, most two-particle interference arrangements 
use entangled photons because of the difficulty in preparing massive entangled 
states. 

The plan of the paper is as follows. First, in Sect. 2, we briefly review 
the basics of decaying systems. The model of interaction, based on coupled 
oscillators, is presented in Sect. 3. In Sect. 4 we derive the states of minimum 
uncertainty, which reach the best compromise between momentum and position 
initial uncertainties. Our analysis is restricted to these states. The evolution 
of the interacting states is evaluated in Sect. 5. Section 6 deals with the 
comparison between perfect correlations and alignment for interacting and free 
systems. In Sect. 7 we evaluate the entanglement present in the system using 
the Schmidt number. Finally, in the Conclusions we discuss the principal results 
of the paper. 

2 Basics of decaying systems 

We present in this Section a brief review of the problem. The classical variables 
describing a decaying process are x + y and p x + P y , with x and y the positions 
of the particles and p x and p y their momenta. In classical theory the law of 
momentum conservation implies that if initially the total momentum is p x + 
p y = 0, it will remain so during the subsequent free evolution of the complete 
system; as p x = p y the particles will always be found in opposite directions. 
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In quantum theory x + y and p x + p y are replaced by the operators x + y 
and p x + p y . The uncertainty relations for these operators are 

A(xi+yi)A(p Xi +p yi )>h (1) 

with i = 1,2,3 the three components of each operator and Axi =< (zj— < 
£i >) 2 > 1 / 2 , • • •. Because of these relations it is not possible to prepare the 
system in a state with perfect correlations in momentum, as it was the case in 
classical theory (A(xi + iji) would have an unbounded value). The most one 
can expect in realistic conditions is to have a sharp distribution centered around 
Px + Py = 0. For instance, in Ref. [2] the wavefunction of the system is taken as 
tp(x, y, t) = J J /(p x , p y ) exp (i(p x • x + p y • y Et)/h)d 3 p x d 3 p y with E the 
energy of the system and / the momenta distribution, sharply centered around 

Px + Py = 0. 

As a consequence of this uncertainty in the initial momenta and the inherent 
spreading of the wavefunction the correlations present in the classical case can 
be lost in many situations in the quantum realm. We discuss this point at extent 
in Sect. 6. 

3 Interaction: coupled oscillators 

Up to our knowledge in all the studies of decaying systems presented so far in the 
literature the subsequent evolution of the system is assumed to be free. However, 
in realistic conditions, in most of the cases there is some type of interaction 
between the products of the process. For instance, if a neutral particle decays 
into two particles with opposite charges there is an electromagnetic interaction 
which, at least at the initial times, is not negligible. The exact description of 
the system taking into account the electromagnetic interaction is, however, too 
complex. 

We propose here a different approach. Instead of studying particular forms 
of realistic interactions we shall analyze an interaction that can be treated an- 
alytically. This approach can be considered as a first approximation to the 
problem. As signaled in the introduction harmonic oscillators are extensively 
used in classical and quantum theory to study complex systems not analytically 
solvable. Then we model the inter-particle interaction via the coupled oscillator 
potential. 

There are other more physical reasons for this choice. Imagine we try to 
introduce some external interaction, for instance through the presence of some 
external fields, aimed to preserve the initial correlations. Examples can be 
found in the literature of the use of external fields to manipulate the form 
of the wavefunction. For instance, an imaginary potential can compensate the 
spreading of the wavefunction resulting in a non-spreading wavepacket [7] . Good 
correlations are obtained when x mea s=a — y mea , i.e., when the measured positions 
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of the particles are close to the perfect (classical) correlation x„ lea = — y. 
Then if we introduce a potential of the type 



nxj)^(x + y) 2 



(2) 



we have an interaction that acts when the particles tend to deviate from the good 
correlations. The strength of the interaction is proportional to the magnitude 
of the deviation, representing the constant k the value of this proportionality. 
This is the ideal form of an interaction aimed to maintain good correlations. It 
can be seen as the ideal limit towards which any realistic interaction with this 
aim must tend. Other acceptable choices should have the form K(x + y) 2 ™, with 
n any positive integer number. However, we shall restrict our considerations to 
the case n = 1. 

The form of Eq. (2) is that of the interaction term of two coupled oscillators. 
Classically, the Hamiltonian of two coupled oscillators is: 



The quantum Hamiltonian is obtained from the classical one by the usual quan- 
tization procedure. Our case corresponds to the choice C\ = C*2 = k and 
C\2 — 2k. As signaled in the introduction coupled oscillators have been used to 
simulate other physical systems [4, 5, 6]. 

The choice of a coupled oscillator to simulate the interaction of the decaying 
particles could be criticized because oscillators are usually associated with states 
of systems that remain bound, whereas decaying systems become free. However, 
it must be noted that the oscillator-type force only acts on the center of mass 
coordinate x + y, remaining the relative position coordinate x y free. This 
way we can have a coupled oscillator interaction and, at the same time, free 
evolution. The particles can become well separated as can easily be tested from 
the solutions of the problem (see comment after Eq. (20)). 

4 States of minimum uncertainty 

As remarked in Refs. [2] and [3] and previous sections the initial uncertainty 
in position and momentum is one of the causes of alignment lost in decaying 
systems. This uncertainty can be minimized choosing the states of minimum 
uncertainty of the relevant operators of the problem, x + y and p x + p y . These 
states can be easily calculated using well-known techniques [8] : given two oper- 
ators Q + and Q_ they are determined by the equations, 



H = 




+ dxj + C 2 *l + Ci 2 x • y 



(3) 



(Q+ - Q+M >= -i<Q- - Q-M > 



(4) 



with 



Q± =< Q± >V> i e = ~ 



<#[Q + ,Q_]|V>> 
2(A V ,Q_) 2 



(5) 
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where [Q+, Q-] is the commutator of Q+ and Q- and A^Q- =< (Q-— < Qi > 
) 2 > x / 2 is evaluated in the state tjj. In the position representation the operators 
can be written as (Q+)j = Xj + yj = xj + yj and (Q-)j — {Px)j + (Py)j = 
—ihd/dxj — ihd/dyj, with j = 1, 2, 3. Then Eqs. (4) and (5) become: 

( Xj + yj - {Q + )i - iiQ-)^ = - ej h (j^- + (6) 



and 

< VI -2h\i> > 



(7) 



where, and from now on, we assume by simplicity e± = €2 = £3 = £■ 

In the particular case Q + = Q- = 0, it is immediate to verify by direct 
substitution that the solution of the above equation is 

V>(x,y) ~exp(-(x + y) 2 /4e?i) (8) 

We note that this is an entangled wavefunction, which reflects a good initial 
correlation for position measurements (x mea — y mea ) with a dispersion of the 
order 2Ti 1/2 e 1 ' 2 . 



5 Evolution of the system 

The evolution of the state is ruled by Schrodinger's equation (we assume by 
simplicity the mass of the two particles to be the same): 

in— ^(x, y, t) = - — (V 2 + V 2 )V(x, y, t) + k(x + y) 2 V(x, y, t) (9) 

In order to solve this equation it is convenient to introduce the center of mass 
and relative position coordinates, defined by the relations: 

X=i(x + y) ; Y = x-y (10) 
The Schrodinger equation becomes: 

ih^(X, Y,t) = — V 2 , + ^ V 2 ,) V(X, Y, t) + ^X 2 ^(X, Y, t) (11) 

where M = 2m is the total mass, /x = m/2 is the reduced mass and ui 2 = 8k/M 
is a habitual way of expressing the coefficient of the interaction potential. 

The solution of this equation at time t can be obtained by integration of the 
initial wavefunction using the kernel or propagator of the system [9] : 

V(X, Y, t) = f d 3 X Q f d 3 Y Q K(X, Y, t; X , Y , i )^ (X , Y Q , t a ) (12) 
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with ^ (X , Y , t ) — N exp (— X 2 /a 2 ) the initial wavefunction, given by 
the minimum uncertainty state studied in Sect. 4. N a represents the normaliza- 
tion of this initial wavefunction (see later) . All the integrations are carried out 
between the limits — oo and oo. As Eq. (11) can be separated for variables X 
and Y, the kernel is the product of the kernels corresponding to the two spatial 
variables [9]: 

K(X, Y, t; Xo, Y„, t ) = K X (X, t; X„, t )K^(Y, t; Y Q , t Q ) (13) 

Now, we have two well-known kernels. From now on we take to = by simplicity. 
For variable X it corresponds to the kernel of the harmonic oscillator [9]: 

Kx(X,t;X ,t ) = ( - .f^ V 7 exp ( %M " ((X 2 + X 2 )cos^ - 2X ■ X Q ) 
\ 2mn sin tot J \2nsmcjt 

(14) 

On the other hand, for variable Y we have a free-particle evolution equation 
with kernel [9]: 

^ y > * y - *•> = ( d*) 3/2 cxp (M (y - Y ° )2 ) (i5) 

After a simple integration using the well-known relation J dzexp(az 2 + j3z) = 
{-it/ - a) 1 / 2 cxp(-/3 2 /4a) valid for Re(a) < 0, we obtain 

^(X,t)=JVexp(i/(X,t)) ( . M ™ -) 1 exp(-a(t)X 2 ) 

\2iTiha sin Lot + Mttlo cos Lot I 



(16) 



where 



M 3 L0 3 I COS LOt \ . MLO COS LOt \ , r2 



and 



/PM) = l^^nT^ U + ^^J + X 

M 2 w 2 a 2 , x 

a(t) = — 5 = 18) 

4^ 2 sin 2 wt + M 2 a 4 w 2 cos 2 tot 

N is the normalization factor. On the other hand, / includes all the real func- 
tions that appear in the form exp(if) in the wavefunction. They are of no 
interest because when calculating probabilities (the magnitudes of interest in 
next section) they only contribute as a constant term. Note that the wavefunc- 
tion does not depend on the Y variable. The initial wavefunction ip a does not 
depend on Y Q . Then as is a free propagator cannot generate a dependence 
on Y. In physical terms, if the initial state does not depend on the relative coor- 
dinates, the interaction term that is only function of the center of mass variables 
cannot introduce that type of dependence during the subsequent evolution. 
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The normalization factor is usually evaluated from the condition for the 
probability / d 3 x J d 3 y|^(x, y, t)\ 2 — 1. However, wavefunction (16) cannot be 
normalized in an absolute sense. Effectively, a simple calculation gives 

J d 3 x J d 3 y|^(x,y,t)| 2 = |Af |M*)| 2 J (19) 

where /j,(t) is the term in Eq. (16) going as the 3/2 power of the expression 
between parentheses. 

Expression (19) is unbounded because the integration is between — oo and oo. 
Then we must work with relative probability densities or probability densities 
per unit volume, which are given by J d 3 x J d 3 y\tp(x, y, t)\ 2 / J d 3 x. It is simple 
to see that this relative probability is normalized to unity. The normalized (in 
this relative sense) wavefunction is: 



3/4 

"J 



V>(X,i)=(^V ox],(-n(/)XN (20) 



Note that there is an additional pure exponential factor exp(i^? M (i)), where </? M is 
the phase of fi(t) (fi(t) — \fi(t)\ exp(iy> M (t))), but it is irrelevant for probabilities 
and can be included in /(X, t). 

The above equation shows that, although we have used an oscillator- type 
interaction, the particles can become well separated. The solution only depends 
on X. The variable Y is not constrained by Eq. (20) and can reach arbitrarily 
large values. 

Now, we consider the decaying system in free evolution. Its mathematical 
form can be derived from the equations in interaction using the relation 

sinwi 

lim = t 21 

First, we note that from the term exp(z/) it is simple to see that the wavefunc- 
tion at times t > has no longer the form of a minimum uncertainty wavepacket 
(as it is well-known [8] , the same behavior occurs for one-particle minimum un- 
certainty wavepackets) . The normalized wavefunction (in the relative sense) is 
given by Eq. 

(20) with a(t) replaced by ap(t), given by 

aF ^ = TTTw (22) 

6 Correlations 

We analyze in this section if the introduction of the interaction improves the 
behavior of the correlations. Two types of measures for the correlations will be 
used. One is the probability of detection with perfect (classical) correlation, the 
other is the alignment at large times. We study them separately. 
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6.1 Perfect correlations 



Two particles have perfect (classical) correlations when the positions of the 
particles found in a measurement process are opposite: x mea = — y mea or 

^med — 0- 

Let us first consider the free evolution case. From the expression for the 
relative probability density for X = we see that the probability of perfect 
correlation goes as a 3 J 2 (t), a decreasing function of time. For very large times, 
ap(t) — > and the probability of perfect correlation becomes negligible. 

We consider now the interacting case, where a(t) shows an oscillatory be- 
havior. It varies between a_ = 1/a 2 and a + = (Mu>a/2h) 2 . Depending on 
the strength of the interaction coupling to we have three different scenarios: 
when cu > 2h/Ma 2 we have a + > a_, whereas for uj < 2h/Ma 2 the relation is 
a + < a_. In the limiting case u> — 2h/Ma 2 we obtain a + = a_. The probabil- 
ity of perfect correlation goes as a 3 ^ 2 (t) and oscillates between a I 2 and a 3 / 2 . 
For instance, for uj > 2h/Ma 2 the probability of perfect correlation increases 
in the interval t € (0, 7r/2) reaching its maximum value at t = tt/2, starting 
a stage of decreasing behavior at this point. The perfect correlations are not 
completely lost even in the limit of very large times, as it was the case for free 
evolution. The analysis for u) < 2h/Ma 2 follows similar lines. The limiting 
situation occurs when uj = 2h/Ma 2 . Now a is constant and the probability of 
perfect correlation detection does not change with time. We conclude that the 
presence of the interaction improves the behavior of the perfect correlations, 
which become an oscillating or constant function instead of the decreasing one 
associated with free evolution. 



6.2 Alignment 

Alignment is a weaker measure of correlations than perfect correlations. We say 
that two particles are aligned when they are detected in, approximately, opposite 
directions. A quantitative criterion for alignment can be introduced following 
Ref. [3] . The angle characterizing angular deviation (perfect alignment is given 
by 7r) can be expressed as: 

TD(t) AX, , s 

tang = — V = — 23) 

R(t) <x,> y ' 

where TD(t) is the transversal deviation and R(t) is the distance both particles 
have traveled. In statistical terms TD can be taken of the order AXi, i. e., the 
variance of any of the components of the center of mass position (i = 1, 2 or 3, 
the problem is isotropic). On the other hand, R is of the order < x% >, the 
expectation value of the position of any of the particles. 

We evaluate now these two variables. We use the Heisenberg picture, where 
the evolution of any operator A (we omit the operator symbol to simplify the 
notation) is ruled by the equation ihdA/dt = [A,H] + ihdA/dt, with H the 
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Hamiltonian of the system. For the position and momentum operators these 
relations become 

dx d y d Px / . x dp y 

m dF^ Px;m dF = p - ^r = - 2K(x+y) = ^r (24) 

Using the center of mass coordinate X and the total momentum P = p x + p y , 
the equations can be rewritten as 

dX 1 dP n v 
m— = -P ; — = -8kX (25) 
dt 2 dt y ' 

The solution of these equations is 

p 

X = X cos Sit + — ^- sin fit (26) 

and 

P = -2mf2X sinQi + P cosQi (27) 

with SI 2 = 4/t/m (note w 2 = 20 2 ) and X Q and P Q the initial values of both 
operators. 

Assuming the expectation values of both initial conditions to be null, < 
X Q >= and < P Q >= 0, the variance of X, can be expressed in the simple 
form: 



(AA^(t)) 2 = (AA M ) 2 cos 2 fit + (AP M ) 2 |^+ < X m P m > ^cobA* 

(28) 

Finally, we evaluate < Xi(t) >. First, we note that from dp x Jdt — — 2/tX, we 
have 

P Xi (t) = P Xi (0) - sin Sit + ^ (cos fit - 1) (29) 

Taking mean values and remembering our assumption < X >= and < 
P >= we have < p Xi (t) >—< p Xi (0) >• Now, we use the relation mdxi/dt = 
p Xi that gives (assuming < Xi{o) >= 0), m < Xi(t) >=< p Xi > t, or 

< Xi(t) >=< p Xi (0) > t/m (30) 
Using again the relation (21), we have for the free evolution 

(AAf {t)f = (AA OT ) 2 + (APo,) 2 ^ + < X oiPoi > (31) 

On the other hand, Eq. (30) is valid for both free and interacting evolutions. 

Now, we can evaluate the angular deviation. We shall concentrate on the 
large times behavior. For the free evolution of the minimum uncertainty state 
we have (the co superscript refers to t — > oo): 

tan 0oo mAP Q* < A > (32) 
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where we have used the relation (1) and the equation < p Xi (0) >= h/ < X >, 
with < A > the mean wavelength associated with the initial mean momentum. 
Equation (32) shows that at large times there is only alignment if < A ><C AX„i, 
i. e., if the mean initial wavelength of the particles is much smaller than the 
dispersion of the center of mass of the source. 

On the other hand, when the interaction is present the angular deviation at 
large times for states with < p Xi (0) is: 

tan0°° = lim = (33) 

t^oo < p Xi (0) > t 

because AXi is a finite function of time due to its periodic behavior. With 
interaction there is alignment at large times independently of the relation be- 
tween the mean initial wavelength and the dispersion of the center of mass of 
the source. 

We remark that result (32) is only valid for minimum uncertainty states (or 
those obeying the approximate relation 2AP iAX i w h), whereas Eq. 

(33) is valid for a much larger class of states, those with < P oi >=< X oi >= 
and <p Xi (0) 0. 



7 Entanglement 

In the previous section we have analyzed the correlations existent in the system. 
Another way to study the problem is to consider the entanglement present in 
the wavefunction. The correlations are the manifestation of the entanglement 
in measurement processes. A measure of the entanglement degree is given by 
the Schmidt number [10, 11, 12]. It has been used in a series of studies with 
important resemblances with our work [11, 12]. In particular, in [11] it was 
studied the entanglement existent between a photon and the atom that has 
emitted it. 

The Schmidt number [10, 11] is given by 

1 1 

S= T^) = T^[pJ) (34) 

with the reduced density matrices 

Px - Tr y {\4, >< VI) ; Py = Tr x (|V >< VI) (35) 

where Tr z denotes trace with respect to the variables associated with particle 
z and |V > is given by Eq. (20). Note that we have used the notation S for the 
Schmidt number instead the usual K one in order to avoid any confussion with 
the propagator. 
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In order to evaluate the above traces we must introduce the Fourier decom- 
position of the wavefunction 

V(x,y,i) - J J dWq^(p,q,i)e ip - x/? V^ = 

d 3 pd 3 q<?!)(p, q,i) < x|p >< y|q > (36) 



where we have used the usual expression < x|p >= (2ttH) 3 / 2 exp(ip • x/h). 
A simple calculation gives 

/ \ 3 / 4 / 2 \ 
0(p,q,t) = ^(-^) exp U 3 (p-q) (37) 



2a(t)J ^\4h 2 a(t) 
On the other hand, we have 

Px = 7V q (|V >< VI) = J d 3 q< q|V >< V|q > = 

d 3 qd 3 pd 3 p*0(p, q,i)0*(p*,q,t) < x|p >< p*|x > (38) 



where the trace is over the momentum q associated with the variable y and 
we have used the relations < q| < y|Q >< Q*|y > |q >=< y|y > <5 3 (q — 

Q)'V'!Q q). 

The second trace is 

Tr P (f>l) = J d 3 p < p| Px |p >= (39) 

d 3 pd 3 qd 3 p»ci 3 q»0(q, p, £)</>* (q, q*, *)</>(p*, q*, t)0*(p*, P, t) 

where the trace now is over the momentum p associated with the variable x. 
Note that in the original approach [10], S counted the number of effective modes 
in the Schmidt decomposition of \ip >■ Here, we do not consider the Schmidt 
decomposition of the wavefunction but the decomposition on the plane wave 
basis. Now the inverse of Tr p (p x ) counts the number of effective modes in this 
basis. This number is still denoted as the Schmidt number (see Ref. [12] where 
this terminology is also used although there is not a Schmidt decomposition of 
the analyzed system). 

A simple calculation gives 

/ \ 3/2 

Tr p (f>l) = V8(2^) 3 (^J (40) 

This equation shows that TY P (/3 X ) has dimensions of momentum 3 longitude 6 . In 
order to obtain a dimensionless Schmidt's number the usual expression above 
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must be replaced. One dimensional parameter provided by the problem is %, 
which suggests to use (2ttU) 3 /Tr p (p^). This expression still has dimensions of 
longitude^ 3 , reflecting the fact that the wavcfunction is normalized as a relative 
probability density per unit volume. Thus, associated with this wavefunction 
we cannot have entanglement but entanglement per unit volume. Then the 
entanglement in a region of volume V is measured by the following generalization 
of Schmidt's number 

Using this equation we obtain 

S v (t) ~ a {tf' 2 (42) 

In particular, in the free case for large times, t » Ma/2%, we have 

(S v )>>(t)~t- 3 (43) 

Equation (42) shows that in the interacting case entanglement is an oscillatory 
function of time. It can be analyzed following the lines of subsection 6.1. This 
behaviour of entanglement explains the persistence of correlations. As a matter 
of fact, it shows the same dependence on time (a(t) 3 ^ 2 ) that perfect correlations. 
On the other hand, in the free case, the entanglement decreases following a 
law of type t~ 3 . After some time almost all the initial entanglement is lost. 
We conclude that the pictures of the system obtained using correlations and 
entanglement are similar. 

Finally, we want to remark that in Ref. [11] the relation between entangle- 
ment and uncertainty relations was studied. A similar analysis could be carried 
along the same lines in our case. However, it would enlarge too much the paper. 
We plan to do it in future work. We only notice that in that reference it was 
also signaled how the uncertainty relations adopt different forms depending on 
the type of measurement done on the system, single-particle or coincidence one. 
In this context it must be remarked that in this work (in Sect. 2) we have only 
considered single-particle measurements, which are usual to the Heisenberg un- 
certainty relations. In a scheme with conditional measurements these relations 
would be different [11]. 

8 Discussion 

We have analyzed in this paper how the behavior of the fragments originated 
by the decaying of a mother particle is modified when interactions between the 
daughter particles or external fields driving the system are taken into account. 
We have shown that the initial correlations present in the system (such as mea- 
sured by perfect correlations and alignment) are much better conserved when 
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the interaction is taken into account. This property, if corroborated for more 
realistic interactions could provide us with sources of particles with good corre- 
lations in those cases where the free evolution degrades the initial correlations. 
In this context note that our analysis also provides a criterion (the relation in 
Eq. (32) between < A > and AX oi ) to determine when alignment is preserved 
by the free evolution. 

Good correlations could be useful for various purposes, for instance, in in- 
terferometric experiments with two-particle massive systems. One example in 
foundational issues is Popper's argument [13, 3], where the actual experiments 
[14] have been carried out with photons instead of massive particles. When the 
experiments are restricted to the framework of quantum optics we lose the possi- 
bility of studying the dependence of the system on the form of the wavefunction, 
which is specially important for interferometric experiments. 

As signaled before decaying systems are archetypal in physics. In the quan- 
tum realm they are particularly well suited to analyze the influence of uncer- 
tainty and spreading on the evolution of the system and, in particular, on its 
correlations. Then the example presented here has also a pedagogical interest. 
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